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Abstract 

The interplay and competition of magnetic and superconducting 
correlations in the weakly interacting two-dimensional Hubbard Model 
is investigated by means of the functional renormalization group. At 
zero temperature the flow of interactions in one-loop approximation 
evolves into a strong coupling regime at low energy scales, signalling 
the possible onset of spontaneous symmetry breaking. This is further 
analyzed by a mean-field treatment of the strong renormalized inter- 
actions which takes into account magnetic and superconducting or- 
der simultaneously. The effect of strong correlations on single-particle 
properties in the normal phase is studied by calculating the flow of the 
self-energy. 



1 Introduction 

The two-dimensional Hubbard model has been proposed by Anderson [Q 
as a basic model for the electronic degrees of freedom in the copper-oxide 
planes of high-temperature superconductors. In agreement with the generic 
phase diagram of the cuprates, the Hubbard model exhibits antiferromag- 
netic order at half-filling, and is expected to become a d-wave superconduc- 
tor near half- filling in two dimensions |2] • The exchange of antiferromagnetic 
spin fluctuations has been suggested as a plausible mechanism leading to d- 
wave pairing [3J |IJ In this picture the BCS effective interaction Vkk' 
is roughly proportional to the spin correlation function x s (k — k'). Close 
to half- filling, x s (q) has a pronounced maximum at the antiferromagnetic 
wave vector (tt, tt). As a consequence, the gap equation with Vkk' as in- 
put has a solution with d x .2_ y 2-wave symmetry, such that the gap Ak has 
maximal modulus but opposite sign near the points (tt, 0) and (0, tt) in the 
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Brillouin zone, respectively. This intuitive argument has been corroborated 
by a self-consistent perturbative solution of the two-dimensional Hubbard 
model within the so-called fluctuation-exchange approximation 

The spin-fluctuation mechanism for pairing could in principle be spoiled 
by other contributions to the BCS interactions and also by spin density wave 
instabilities. Unfortunately, it is hard to prove the existence of supercon- 
ductivity in the Hubbard model by exact numerical computation j^l Ej , as a 
consequence of finite size and/or temperature limitations. However, the ten- 
dency towards antiferromagnetism and d-wave pairing is captured already 
by the 2D Hubbard model at weak coupling, where one can hope to pro- 
ceed by perturbative calculations. Conventional perturbation theory breaks 
down for densities close to half-filling, since competing infrared divergences 
appear as a consequence of Fermi surface nesting and van Hove singularities 
[HI El IH3 • A controlled and unbiased treatment of these divergencies cannot 
be achieved by standard resummations of Feynman diagrams, but requires a 
renormalization group (RG) analysis which takes into account the particle- 
particle and particle-hole channels on an equal footing. In the following we 
describe progress in this direction achieved by applying the socalled func- 
tional renormalization group (fRG) to the two-dimensional Hubbard model. 

In Sec. 2 we define the two-dimensional Hubbard model and describe the 
relevant infrared singularities in perturbation theory. In Sec. 3 we briefly 
review the Wick ordered version of the functional renormalization group and 
discuss the one-loop flow of effective interactions and susceptibilities for the 
Hubbard model. In Sec. 4 we compute the flow of the electronic self-energy 
and discuss how the strong correlations at low energy scales affect spectral 
properties of single-particle excitations. In Sec. 5 we treat the effective low- 
energy problem obtained from the fRG flow at zero temperature within a 
mean-field approximation and find a range of fillings where magnetic and 
and sizable superconducting order coexist. A short conclusion follows in 



2 Two-dimensional Hubbard Model 

The Hubbard model describes tight-binding electrons with a local repulsion 
U > 0. The Hamiltonian reads 



where c\ a and c- lrT are creation and annihilation operators for fermions with 
spin projection a =1,1 on a lattice site i, and nj CT = cj^cjo-. A hopping 
amplitude —t between nearest neighbors and an amplitude —if between next- 
nearest neighbors on a square lattice leads to the dispersion relation 
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(1) 



ek = — 2i(cos k x + cos k y ) — 4t'(cos k x cos k y ) 



(2) 
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for single-particle states. This dispersion relation has saddle points at k = 
(0,7r) and (vr, 0), which lead to logarithmic van Hove singularities in the 
non- interacting density of states at the energy e v H = 4i'. 

For a chemical potential \i = e v H the Fermi surface contains the van Hove 
points. In this case a perturbative calculation of the two-particle vertex 
function leads to several infrared divergencies already at second order in U, 
that is at one-loop level [HI El El- In particular, the particle-particle channel 
diverges as log 2 for vanishing total momentum ki + k2 , and logarithmically 
for ki + k2 = (it, it). The particle-hole channel diverges logarithmically 
for vanishing momentum transfer; for momentum transfer (tt, tt) it diverges 
logarithmically if t' 7^ and as log 2 in the special case if = 0. Note that 
\i = for if = corresponds to half-filling (n = 1). In this case there are 
also logarithmic divergences for all momentum transfers parallel to (tt, tt) 
or (tt, —tt) due to the strong nesting of the square shaped Fermi surface. 
For fi ^ e v H only the usual logarithmic Cooper singularity at zero total 
momentum in the particle-particle channel remains. However, the additional 
singularities at fj, = e v H lead clearly to largely enhanced contributions for 
small \fi — e v n\, especially if t' is also small. 

Hence, for small \fi — e v n| one has to deal with competing divergencies in 
different channels. This problem calls for a renormalization group treatment. 

3 Renormalization Group Approach 

Early RG studies of the two-dimensional Hubbard model started with simple 
scaling approaches, very shortly after the discovery of high-T c superconduc- 
tivity [HI El E31- These studies focused on dominant scattering processes 
between van Hove points in k-space, for which a small number of running 
couplings could be defined and computed on one-loop level. Spin-density 
and superconducting instabilities were identified from divergencies of the 
corresponding correlation functions. 

A complete treatment of all scattering processes in the Brillouin zone is 
rather complicated since the effective interactions cannot be parametrized 
accurately by a small number of variables, even if irrelevant momentum and 
energy dependences are neglected. The tangential momentum dependence 
of effective interactions along the Fermi surface is strong and important in 
the low-energy limit. Hence, one has to deal with the renormalization of 
coupling functions. This problem is treated most naturally within a Wilso- 
nian momentum shell RG where modes are integrated out successively 
from high energy states down to the Fermi surface [121 ■ This type of RG 
is also the basis for important rigorous work on two-dimensional Fermi sys- 
tems ^Hl- The successive integration of modes can be formulated as an ex- 
act hierarchy of flow equations for effective interactions (one-particle, two- 
particle etc.), which is known as the exact or functional renormalization 
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group |14[I15[ IT6]. Starting with the work of Zanchi and Schulz |17j . several 
groups have computed effective interactions and susceptibilities for the two- 
dimensional Hubbard model, using various versions of the fRG in one-loop 
approximation |18l I19j . Here we focus on results obtained from the Wick- 
ordered version of the fRG Other versions yield qualitatively similar 
results. 



3.1 Wick ordered flow equations 

The renormalization group equations are obtained as follows (for details, see 
Salmhofer |13j and Ref. JH])- The infrared singularities are regularized by 
introducing an infrared cutoff A > into the bare propagator such that con- 
tributions from momenta with |ek — n\ < A are suppressed. All Green func- 
tions of the interacting system will then depend on A, and the true theory is 
recovered only in the limit A — > 0. The RG equations are most conveniently 
obtained from the effective potential V A , which is the generating functional 
for connected Green functions with bare propagators amputated from the 
external legs. Taking a A-derivative one obtains an exact functional flow 
equation for this quantity. Expanding V A on both sides of the flow equation 
in powers of the fermionic fields (i.e. Grassmann variables), and comparing 
coefficients, one obtains the so-called Polchinski equations |15) for the effec- 
tive m-body interactions used by Zanchi and Schulz ^7]. Salmhofer ^3] has 
pointed out that an alternative expansion in terms of Wick ordered mono- 
mials of fermion fields yields flow equations for the corresponding m-body 
interactions V A with a particularly simple structure (see Fig. 1). The flow 
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Figure 1: Diagrammatic representation of the flow equation for in the 
Wick ordered version of the fRG. The line with a slash corresponds to 
dD A /dA, the others to D A ; all possible pairings leaving m incoming and 
m outgoing external legs have to be summed. 



of Vj^ is given as a bilinear form of other n-body interactions (at the same 
scale A), which are connected by lines corresponding to the propagator 
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(3) 



where £k = £k — A*> and one line corresponding to dD A (k)/dA. For small A, 
the momentum integrals on the right hand side of the flow equation are thus 
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Figure 2: Fermi surfaces (solid lines) and support of the propagator D A 
(shaded regions) in momentum space; left: t 1 = 0, right: t' > 0. The Fermi 
surface on the right intersects the socalled umklapp surface (dashed line). 

restricted to momenta close to the Fermi surface (see Fig. 2). With the ini- 
tial condition V A ° = bare interaction, where Ao = max |£k|, the above flow 
equations determine the exact flow of the effective interactions as A sweeps 
over the entire energy range from the band edges down to the Fermi surface. 
The effective low-energy theory can thus be computed directly from the mi- 
croscopic theory without introducing any ad hoc parameters. The structure 
of these flow equations is very convenient for a power counting analysis to 
arbitrary loop order |13| . and also for a concrete numerical solution. 

3.2 One-loop flow 

To detect instabilities of the system in the weak-coupling limit, it is suffi- 
cient to truncate the infinite hierarchy of flow equations described by Fig. 
1 at one-loop level, and neglect all components of the effective interaction 
except the two-particle interaction V^- The effective two-particle interac- 
tion V^ A then reduces to the one-particle irreducible two-particle vertex T A , 
and its flow is determined by T A itself (no other m-body interactions enter). 
Putting arrows on the lines to distinguish creation and annihilation oper- 
ators one thus obtains the flow equation for T A shown graphically in Fig. 
3. Flow equations for susceptibilities are obtained by considering the exact 
fRG equations in the presence of suitable external fields, which leads to an 
additional one-body term in the bare interaction, and expanding everything 
in powers of the external fields to sufficiently high order ^Sj. On one-loop 
level one obtains the flow equations shown in Fig. 4. The flow of a suscepti- 
bility x A is determined by the corresponding response vertex the flow of 
which is in turn given by T A and i? A itself. The initial conditions are given 



5 




1 PH 2 1 PH' 2 



Figure 3: Flow equation for the two-particle vertex T in one-loop approx- 
imation with the particle-particle channel (PP) and the two particle-hole 
channels (PH and PH'). 




Figure 4: Flow equations for (a) the susceptibilities x an d (b) the response 
vertices in one-loop approximation. 

by X A ° = f° r the susceptibilities and by the bare response vertices for R A . 
For pairing susceptibilities only the particle-particle channel contributes to 
the propagator pair in Fig. 4, for charge and spin density susceptibilities 
only the particle-hole channel. 

It is clearly impossible to solve the flow equations with the full energy 
and momentum dependence of the two-particle vertex, since T A has three 
independent energy and momentum variables. The problem can however be 
simplified by ignoring dependences which are irrelevant (in the RG sense) 
in the low energy limit, namely the energy dependence and the momentum 
dependence normal to the Fermi surface. Hence, we compute the flow of 
the two-particle vertex at zero energy and with at least three momenta on 
the Fermi surface (the fourth being determined by momentum conservation) . 
On the right hand side of the flow equation we approximate the vertex by its 
zero energy value with three momenta projected on the Fermi surface (if not 
already there), as indicated in Fig. 5. This projection procedure is exact for 
the bare Hubbard vertex, and asymptotically exact in the low-energy regime, 
since only irrelevant dependences are neglected. The remaining tangential 
momentum dependence is discretized. The momentum-dependence of the 
two-particle vertex is thus approximated by a step function which is constant 
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Figure 5: Projection of momenta on the Fermi surface; discretization and 
labelling of angle variables. 



on "patches" (sectors) in the Brillouin zone. 

The flow of the two-particle vertex has been computed for many differ- 
ent model parameters t' and U (t just fixes the absolute energy scale) and 
densities close to half-filling. At T = the vertex develops a strong momen- 
tum dependence for small A with divergencies for several momenta at some 
critical scale A c > 0, which vanishes exponentially for U — > 0. To see which 
physical instability is associated with the diverging two-particle vertex, var- 
ious susceptibilities have been computed, in particular: commensurate and 
incommensurate spin susceptibilities Xs(q) with q = (7r,7r), q = (vr — 5, n) 
and q = (1 — 5)(tt,tt), where 5 is a function of density j^Jj, the commen- 
surate charge susceptibility xc(( 7r ; n)), and singlet pair susceptibilities with 
form factors 



d(k) = < 



1 (s-wave) 

^ (cos k x + cos k y ) (extended s-wave) 

■^(cos k x — cos k y ) (d-wave d x i_ y i) 

sin k x sin k y (d-wave d xy ). 



(4) 



Some of these susceptibilities diverge together with the two-particle vertex 
at the scale A c . Depending on the choice of U, t' and \x the strongest diver- 
gence is found for the commensurate or incommensurate spin susceptibility 
or for the pair susceptibility with d x 2_ y 2 symmetry. In Fig. 6 we show a 
typical result for the flow of the two-particle vertex and susceptibilities as 
a function of A. Only the singlet part (from a spin singlet-triplet decompo- 
sition) of the vertex is plotted, for various choices of the three independent 
external momenta. The triplet part flows generally more weakly than the 
singlet part. Note the threshold at A = 2|/x| below which the amplitudes for 
various scattering processes, especially umklapp scattering, renormalize only 
very slowly. The flow of the antiferromagnetic spin susceptibility is cut off 
at the same scale. The infinite slope singularity in some of the flow curves at 
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Figure 6: (a) The flow of the singlet vertex as a function of A for several 
choices of the momenta kpi, \<lf2 an d which are labelled according to 
the numbers in Fig. 5. The model parameters are U = t and t' = 0, the 
chemical potential is \x = — 0.02t; (b) the flow of the ratio of interacting and 
non-interacting susceptibilities, X^/Xo, f° r the same system [from Ref. |18]]. 

scale A = \fj,\ is due to the van Hove singularity being crossed at that scale. 
The pairing susceptibility with d x 2_ y 2 -symmetry is obviously dominant here 
(note the logarithmic scale). Following the flow of the two-particle ver- 
tex and susceptibilities, one can see that those interaction processes which 
enhance the antiferromagnetic spin susceptibility (especially umklapp scat- 
tering) also build up an attractive interaction in the d x i_ y i pairing channel. 
This confirms the spin-fluctuation route to d-wave superconductivity |2j. 
Similar results leading to the same conclusions have been obtained by other 
one-loop calculations using different versions of the fRG jl3 El [20j • 
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4 Two-loop self-energy 



The self-energy £ for the 2D Hubbard model, which encodes the influence 
of correlations on single-particle excitations, has been analyzed within the 
functional RG framework in several recent articles. Interesting dynamical 
contributions to the self-energy set in at second order in the two-particle 
interactions, corresponding to the two-loop diagram shown in Fig. 7. The 



vertices on the right-hand side were obtained from the one-loop flow equa- 
tion. A full two-loop calculation, including two-loop vertex renormalization, 
would be very complex, and has not been done so far. 

Some of the earlier studies focused on the frequency derivative of the 
self-energy on the Fermi surface, which determines the wave function renor- 
malization factor Z = [1 — 9 a; ReS(ki?, a;)!^^]" 1 . In Fermi liquids Z yields 
the weight of the quasi-particle peak in the spectral function. Zanchi [22] 
computed the flow of Z for the 2D Hubbard model with pure nearest neigh- 
bor hopping at half-filling and found that Z is strongly reduced when the 
renormalized interactions increase, the suppression being strongest near the 
van Hove points. Subsequently, Honerkamp and Salmhofer [23] computed 
the flow of the Z-factor at and also away from half-filling. While confirming 
an anisotropic reduction of Z, they observed that this suppression evolves 
much slower than the divergence of dominant interactions. In addition, 
they pointed out that the self-energy does not significantly affect the vertex 
renormalization in the perturbative regime, that is before the interactions 
become very strong. The imaginary part of the self-energy at zero frequency, 
lm£(kp,0), was computed by Honerkamp [21]. In Fermi liquids this quan- 
tity is directly related to the decay rate of quasi-particles. For a hole-doped 
Hubbard model with a concave Fermi surface that almost touches the van 
Hove points (this requires a negative t') he found that Im£(kF,0) is mod- 
erately anisotropic, with larger values for momenta closer to the van Hove 
points. 

Functional RG calculations of the full frequency dependence of the self- 
energy were performed recently by Katanin and Kampf [25] within the one- 
particle irreducible version of the fRG, and by two of us [^Hl within the Wick 
ordered version. The real and imaginary parts of T,{kp,uj) and the resulting 
spectral function A(\<.f,uj) were computed for the 2D Hubbard model with 
nearest and next-to-nearest neighbor hopping at finite temperature, which 
was chosen such that the largest renormalized interactions were strongly 
enhanced, but did not diverge for A — > 0. Marked deviations from Fermi 




Figure 7: Flow equation for the two-loop self energy. 
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liquid behavior were obtained for kp close to a van Hove point, and also 
near other hot spots [27] in the case of filling factors above the van Hove 
singularity: the imaginary part of £(kp,u;) develops a pronounced peak at 
low frequencies, as seen for example in Fig. 8, which leads to a double-peak 
structure in the spectral functions, reminiscent of a pseudogap. This double 
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Figure 8: Imaginary part of the self energy ImS(k_p, u) normalized by U 2 as 
a function of u and on a hot spot, for parameters t'/t = —0.1, T/t = 0.05, 
U/t = 1.8 and density n = 0.92; the result from the functional RG (solid 
blue line) is compared to the result from second order perturbation theory 
(dashed black line) [from Ref. [2*S]]. 

peak gradually transforms into a single peak when moving away from the 
hot spot or van Hove regions, and also upon raising the temperature. In 
the one-particle irreducible scheme the peak in Im£(kp, u) was found to be 
partially split by a dip upon moving away from the van Hove point, which 
leads to a three peak structure in the spectral function No such splitting 
was obtained in the Wick ordered functional RG j^H]- This discrepancy 
can be traced back to differences in the approximate flow equations of both 
schemes, and indicates that higher order corrections become important when 
the renormalized interactions grow too strong. 

One can conclude that self-energy contributions are particularly large 
near hot spots and tend to destroy Fermi liquid behavior, but the final form 
of the spectral function requires an analysis beyond the above perturbative 
truncation of the flow equations. 

5 Flow to strong coupling: RG + Mean-Field 

Within the one-loop truncation the renormalized two-particle interaction T A 
always diverges at T = in one or several momentum channels at a finite 
energy scale A c , even for a small bare interaction U, resulting in a strong 



10 



coupling problem in the low-energy limit. If the two-particle vertex diverges 
only in the Cooper channel, driven by the particle-particle contribution to 
the flow, this strong coupling problem can be treated by exploiting A c as a 
small parameter |28j . The scale A c is exponentially small for a small bare 
interaction. The formation of a superconducting ground state can then be 
described essentially by a BCS theory with renormalized input parameters. 
A method very similar in spirit has been applied very recently to supercon- 
ductivity in quasi-one-dimensional systems 29 . 

Here, we use an analogous approach to treat several order parameters 
simultaneously at mean-field level. Near half filling, for U = 2.5t and 
t' = — 0.2t the one-loop flow produces strong interactions in the Cooper 
channel as well as in the particle-hole channel, the latter being dominated 
by commensurate (tt, it) scattering. We thus concentrate on d-wave super- 
conductivity and antiferromagnetism, which can be treated within a 4 x 4 
Nambu formalism. While the interplay between antiferromagnetic order and 
d-wave superconductivity has been studied earlier within mean-field theories 
using models containing pairing interactions |3L)l 131 j , we emphasize that we 
do not introduce such couplings ad hoc. Rather, pairing interactions arise 
within the purely repulsive Hubbard model via the fRG treatment. 

In an extended RG framework, spontaneous symmetry breaking can be 
handled by adding an infinitesimal symmetry breaking term at the beginning 
of the flow, which is then promoted to a finite order parameter at the scale 
A c |32j . However, this approach would be extremely involved for the case 
of competing order parameters we are interested in here. Instead, we stop 
the one-loop flow at a scale Ai that is small compared to the band width 
but still safely above the scale A c where the two-particle vertex diverges. 
At this point the vertex has developed already a pronounced momentum 
dependence, reflecting in particular magnetic and superconducting correla- 
tions. The integration over the remaining modes, below Ai, is treated in a 
mean-field approximation allowing antiferromagnetic and superconducting 
order. The mean-field theory is defined on a restricted momentum region 
near the Fermi surface, with |ek — (j\ < Ai, and the effective interactions 
entering the mean-field equations are extracted from Details on this 

combined RG+MF approach will be presented elsewhere [HS]. We now dis- 
cuss some results for the order parameters. 

In Fig. 9 we show the amplitudes of the d-wave superconducting gap 
(SCG) and the s-wave antiferromagnetic order parameter (AFG) as a func- 
tion of (electron) density. The latter is computed simultaneously with the 
order parameters, and thus differs from the initial density of the bare sys- 
tem. For n = 1, i.e. at half filling, the AFG is largest, while the SCG is 
zero. In this case, the system is fully gapped. With decreasing filling the 
AFG decreases initially roughly linearly, while the SCG remains numerically 
zero. For electron densities below one, holes appear first in pockets around 
(7r/2,7r/2), which define a surface of low-energy excitations of the non-half 
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Figure 9: d-wave superconducting order parameter and antiferromagnetic 
order parameter (in units of t) as a function of density for U/t = 2.5 and 
t'/t = —0.2. Filled colored symbols represent the results obtained for the 
combined theory with two order parameters, while in the calculation leading 
to the open symbols only one order parameter, either antiferromagnetic or 
superconducting, was allowed in the mean-field equations. 



filled system which increases with larger hole doping. In principle, this 
residual Fermi surface is always unstable against superconducitivity, due to 
attraction in the d-wave Cooper channel. However, very close to half-filling, 
where the pockets are small, the SCG is tiny, since the d-wave attraction 
is very small near the Brillouin zone diagonal. When the hole pockets are 
large enough, roughly at n « 0.88, a sizeable SCG develops, coexisting with 
a finite AFG. The numerical analysis is rather involved in the coexistence 
region, and the precise behavior of the order parameters in this regime re- 
mains to be clarified. What is clear is that there is a range of densities for 
which both order parameters are sizable. The coexistence disappears below 
n as 0.85. Below this density the AFG truly vanishes, and the SCG is finite 
and decreases with decreasing filling. The respective results for both order 
parameters when the other one is set to zero are also shown. When the 
AFG is set to zero, the SCG persist even at half filling. When the SCG is 
set to zero, the AFG is enhanced in the coexistence region. In both cases, 
a finite value for one order parameter leads to a suppression of the other. 
While these results have been obtained for a weak on-site repulsion, they are 
in line with the behaviour at stronger coupling as found within variational 
cluster perturbation theory [SI]. There, however, due to finite size effects 
the possibilty of coexistence on the hole doped side could not be conclusively 
answered. 

Fig. 10 shows the angular dependence of the SCG and the AFG for three 
selected densities, along with the corresponding hole pockets. Here, (ft is the 
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Figure 10: Top: d-wave superconducting and antiferromagnetic order pa- 
rameters at three different values of the chemical potential. Bottom: Fermi 
surfaces in the magnetic Brillouin zone (solid blue lines); in the antiferro- 
magnetic state close to half-filling the Fermi surface forms a hole pocket 
around (tt/2, it/2). The corresponding bare Fermi surfaces, backfolded with 
respect to Umklapp surface, are shown as broken lines. The straight lines 
indicate the patching scheme. 

symmetry and is slightly anisotropic. The SCG is practically zero. The oval 
pocket which opens around (tt/2,tt/2) is sizeable, but does not extend far 
enough to allow for a sizeable d-wave SCG. For n = 0.8717 we observe a 
coexistence of both order parameters. While the AFG is reduced in size 
compared to the case n = 0.9453, its shape remains essentially the same. 
The pocket extends further away from the zone diagonal, and together with 
the smaller AFG this allows for a substantial d-wave SCG. In this situation 
the low-energy excitations are gapped due to the AFG near the points (n, 0) 
and (0, 7r), and gapped due to the SCG along the "Fermi pocket" , with nodes 
along the Brillouin zone diagonals. For an even smaller density of n = 0.8228 
the AFG vanishes and only a d-wave SCG remains, which extends over the 
whole Fermi surface, except at the nodal points. 

While we have restricted our analysis to superconductivity and antiferro- 
magnetism, there are other instabilities which may arise, such as ferromag- 
netism (at moderate \t'/t\ in the Hubbard model) [SSI dHl GUI or a d-wave 
Pomeranchuk instability of the Fermi surface |38l l3*§] . The treatment pre- 
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sented here neglects order parameter fluctuations, which are particularly 
important in low dimensions. Since we work at T = we expect the above 
results to be qualitatively stable. In a more sophisticated approach order 
parameter fluctuations can be treated most conveniently by introducing ap- 
propriate bosonic fields, as discussed recently for antiferromagnetic order in 
the half-filled Hubbard model at finite temperature |4()j . 



6 Conclusion 

In summary, the functional renormalization group captures many aspects of 
the complex interplay of magnetic and superconducting correlations in the 
two-dimensional Hubbard model. In particular, it captures the generation of 
attractive d-wave Cooper interactions via antiferromagnetic spin correlations 
already in a one-loop approximation of the vertex flow. Strong correlation 
effects in the single particle excitations are obtained from the two-loop flow 
of the self-energy. Magnetic fluctuations suppress spectral weight at hot 
spots, leading to a pseudogap shape of the spectral function at these points. 
The strong-coupling problem posed by the divergence of the flow at T = is 
treated approximately by combining the fRG with a mean-field treatment at 
low energies. Superconducting and antiferromagnetic instabilities are found 
to compete, with a small range of densities for which both orders can coexist 
with a sizable order parameter for each. 

A more complete analysis of the two-dimensional Hubbard model by 
renormalization group methods, including also order parameter fluctuations 
at low energy scales, can be expected to yield further important clues for a 
better understanding of magnetic and superconducting correlations in high 
temperature superconductors. 
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